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During recent years, the energy changes of high energy partons traversing quark-gluon 
plasma(QGP) have been the subject of intensive interest 0. To study the behavior of the 
high energy partons through QGP, one must constitute a kind of response theory for QGP 
to external current. The color permittivity, characteristic quantity of QGP media, must be 
known. In the frame of kinetic theory, the Abelian-like color permittivity of QGP under the 
linear approximation has been given out long ago @. Just as Ref. || pointed out, the self- 
coupling terms of color field have not been included in the calculation of relevant quantities. 
Ref. [[| gave out the nonlinear color permittivity of QGP beyond the linear approximation, 
in which the SU(N C ) color algebra has been ensured in the iteration process and, to some 
extent, the non-Abelian characteristic of QGP has been reflected. However, there still exists 
the problem pointed out by Ref. 

In this letter, the color permittivity of QGP including the self-coupling contribution of 
the color field is calculated with the kinetic theory. Recently, a few works have been done by 
separating the distribution functions and color field into regular and fluctuation parts, with 
which the semi-classical kinetic equations are converted to the forms which can be solved 
conveniently The quantities in the kinetic equations of QGP, such as the quark or 

antiquark distribution function f q (x,p) or fg(x,p), gluon distribution function G(x,p) , field 
A(x) and induced color current j in QGP, can be decomposed as 

h = f* + fl A, = A« + Al, j=f+j R , (1) 

where the index R represents the regular parts of corresponding quantities, i.e., f R = (f), 
with index T the stochastic fluctuation parts. (■ ■ •) represents taking average over a statistical 
ensemble. Still further, the distribution functions /J,..., can be expanded as the series of 
small fluctuation A T , for example 

oo 

fJ = T,fI H - (2) 

n=l 

With this kind of decomposition, one can obtain the kinetic equations for the fluctuation 
parts of quark, antiquark and gluon distribution functions from the semi-classical kinetic 
equations of QGP, respectively @. If considering only the fluctuations around the equilibrium 
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state, fq{fq)-, G R can be chosen as Fermi-Dirac and Bose-Einstein equilibrium distribution 
functions, respectively. One can set A R = and need only to consider the evolution equation 
of each fluctuation quantity. For example, the fluctuation part of quark distribution function 
satisfies 

ig{K> f q \ - K, f q ])) - \9P»{FU' ^> + \i9 2 P»{[Al Al\ - {[Al A T V \), g} 

-\^{f^ g) - g») + ^(W d f) - mi ah % }))• 

(3) 

The kinetic equations of the fluctuation parts of the antiquark and gluon distribution func- 
tions are similar to Eq.([5]) except the opposite signs of the terms related to {...,...} for an- 
tiquarks; /, A and F^ u are replaced by G, A and T^ v for gluons, respectively. = F® u t a , 
T^y = F® u T a , with A = A a t a and A = A a T a . The t a and T a are the corresponding gener- 
ators of SU(N C ) in fundamental and adjoint representation, respectively. In the following 
discussion we will only give out explicitly the relevant equations for /J, but not for fT and 
G T unless specialized otherwise. In Eq.(^), the index L denotes the linear term of F^ with 
respect to AF. 

The fluctuation part A T of color field A obeys the following Yang-Mills equations 

dX" V = ~f v + W[A r ", A Tv \ - ([A T », A Tu ])) + 
+i9([A T „ if H - ([Al F L r H)) + 9 2 ([Al, [A TfJj , A Tu \] - (LAj, [A Tfl , A Tv ]])), (4) 

and j Tu can be expanded according to Eq.(0), with 

fM» = g? J j0-5P U i trta (fI (n) ~ f?^) + Tr{T a G T ^)l (5) 

being the n-th order color current. 

It is convenient to work in the temporal axis gauge, i.e., A = and in the momentum 
space. The index T of field will be omitted as well, if without confusion. By substituting the 
n-th order distribution function f^ n \k,p), ■ ■ ■ obtained from the n-th order kinetic equations 
into Eq.(|5p, one can obtain j T( - n \k). If expanding j T (k) up only to the first order current 
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j T ^(k), one can recover the Abelian-like permittivity El- It can be confirmed that the 
second order color current does not contribute to the nonlinear permittivity 00. So, to 
obtain the nonlinear permittivity including the self-interaction of gauge field, one needs to 
analyze the third order equations directly. 
The following shorthands 



E =[5(k-h- k^dhdh-^—, M eq = ~ (4V) + /f (p )) + N c G^(p°). 

k=ki+k2 v / 

=Js(k-k 1 -k 2 - h)dhdk 2 dhj^-, X al3 (k,p) = (p • k)g af3 -p a k^. (6) 
are introduced for notational simplicity. The third order kinetic equation is as following 

-i^k^\k,p) = 

*gp l E [([A(k 1 ),f^ 2 \k 2 ,p)}-([A l (k 1 )j^ 2 \k 2 ,p)]))+ 

k=ki+k2 

+\w E x iX (h, P )^: ({A^j^ik^p)} - ({Mh),ff 2) (h2,p)})) + 

+ E V 1 ^ ({[A(fci),^(fe)],/ 9 T(1) (fe,p} - ({[^(^),^(A: 2 )],/ g T(1) (^,p)}) 



(7) 



Considering the soft excitation carrying momentum k ~ gT and the momentum of particles 
p ~ T in QGP IB, by substituting the third order distribution functions (0) into Eq.fl5|), 
and keeping only the leading order in g, one may identify after some lengthy calculation 

f^ a \k) « 
where 

,(J)ii« 4 V- rabcrcde f ^P p l p>p k p l 1 UJ 2 d° p Ka 



^(I)ijkl _ 4 V"^ rabcrcde f 



k=ki+k?+k 3 



(2tt) 3 p° p ■ k + ip°ep ■ (ki + k 2 ) + ip°ep ■ k 2 + ip°e 

(9) 



Now turn to the calculation of the non-Abelian permittivity from the mean field equation. 
By expanding j T (k) up to third order, substituting j T ^ % {k) and j T ^ l (k) leading order in g 
into Eq.(f|) in temporal axis gauge and in momentum space, one can identify 



-gH a 



E \9 

k=ki+k 2 +k 3 V 



d?p 



abc cede 



fabCf 



X 



(2ir) 3 p° p ■ k + ip°e p ■ (k\ + k 2 ) + ip°e 

p • k p • ki p ■ k 2 p • k 3 u 2 dlN eq 
~K K y ~K 2 K 3 p-k 2 + ip°e 

x (A^kJA^A'fa) - A b {k 3 )(A d (ki)A e {k 2 )) - (A b (k 3 )A d (h)A e (k2)) 
k, • k 2 k • k 3fabcfcde \ Ab{h)A d {ki)A e {k2) _ {A ^ h)A ^ kl)A e {k2)) " 



(10) 



K X K 2 KK 3 

The first term on the right-hand side of above equation with the factor in big parentheses 
corresponds to the interaction between the particles and the secondary waves resulting from 
the nonlinear interactions of the eigenwaves in QGP. It should be stressed that the second 
term on the right-hand side of Eq.fllOD with the factor in square brackets is the plasmon 
self-interaction term in QGP, which is characteristic of the QCD plasma and different from 
the QED. It is well known that the solution of the Yang-Mills equation in QCD vacuum is 
very difficult because of the self-coupling term. Just as Ref. || pointed out that, in previous 
works, this term reflecting the non-Abelian characteristics has been discarded. However, 
in hot dense medium where the phase of the excitation of field is random, one can take 
average over phase in the process of calculation @ and make it possible to solve the mean 
field equation (|) or (|H}). 

Multiplying Eq.© with A 9 (k'), then taking average with respect to the random phase, 
one can identify % 

(e L 5 ad + e N ad L )(dg)=Q ) (11) 

with the notation (dg) = (A d A 9 )k representing the correlation strength of color field. The 
nonlinear permittivity defined in Eq.(|ll|) depends on the correlation strength 



.NL 



= c° -l-p 
'ad — °ad 1 ^adi 

,2 



(12) 



'ad 



'ad 



g~ 



or 



f abc f cfd (bf) kl + 



[kk, 

2 / 



{fab e f ede {be)ki+f ade f efe {fe)kA 



1 



(2tt) 3 p° 
1 



K J p ■ k + ip°e 



X- 



p ■ (k — ki) + ip°e \p ■ k + ip°e 



f cfd (bf) kl + 



f cde (be) kl 



p ■ k\ + ip°e 



dki 



dki 



The ef d comes from the self-interaction of plasmons, while the e^ d represents the contribu- 
tion from the the interactions between the plasma particles and the beating of two timelike 
eigenwaves. The factors (p ■ k + ip°e) , (p • k\ + ip°e) 1 , (p ■ (k — ki) + ip°e) 1 in e^ d have 
different contributions to the integrals. The two former factors have only real contribu- 
tions because the eigenwaves in QGP are always timelike @, while the latter has both real 
and imaginary contributions to e R d because the beating of two timelike eigenwaves can be 
spacelike by noticing the crucial relation 



p ■ (k — k\) + ip°e p-{k — k\) p° 1 p° 

From Eq. (|T2"D , one can conclude that the nonlinear permittivity e^ d is a matrix in color 
space and depends on the correlation (A b A c )k- A rough value estimate can make us see the 
contributions of ef d and e^ d to the color permittivity more clearly. One can easily calculate 
the diagonal elements that is considering only the correlation of the same color by using 

7T 

{bc) k = -(5(u - u k ) + 5(u + u k ))I k 5 bc , (14) 

ur 

where Jk characterizes the total intensity of the fluctuating oscillation with frequency u^. 
Furthermore, if considering the fluctuation of thermal level, i.e., one can take Ik = 4nT. By 
taking the upper limit of K\ to be gT while the lower limit g 2 T, the real part of the diagonal 
elements of e^ d can be easily obtained for the important k = modes of the pure gluon 
plasma occasion 

Re(e^ L ) = 5 ad {e s + Ree R ) w -5 ad (0.9763 + 1.913#). (15) 

Summarizing, the non-Abelian permittivity ([12]) including the self-interaction of gauge 
field is given out for the first time. It relies on the concrete modes a>k and the correlation 
of field as in the electro-magnetic plasma. However, It is truly non-Abelian, because of the 
presences of color indices and the SU(N C ) structure constants. The motivations for us to 
calculate the permittivity including the self-interaction are the following: 1. The gauge in- 
variance is a embarrassing and important problem. Without considering the self-interaction, 
the SU (N c ) gauge symmetry will be violated. The color permittivity reflecting the character- 
istic of QGP must include the influence of self-interaction. 2. The non-Abelian permittivity 
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contributed by the self-interaction will influence the relevant physical quantities of QGP. For 
example, the rough numerical result fll"5|) indicates that the nonlinear eigenfrequency shift is 
changed greatly by this effect H 



-5 



= = 0.282^T. 
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